Ferrotoroidic Moment as a Quantum Geometric Phase 
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We present a geometric characterization of the ferrotoroidic moment r in terms of a set of Abelian 
Berry phases. We also introduce a fundamental complex quantity, 2^^, which provides an alternative 
way to calculate r and its moments and is derived from the tensor T^j^ = 2 ^j^j- This geometric 
framework defines a natural computational approach for density functional and many-body theories. 
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Introduction. The recent discovery of new compounds 
with strong multiferroic couphng revised interest in mul- 
tiferroic materials. An attractive component of this dis- 
covery is the potential for realizing strong magnetoelec- 
tric effects, i.e., magnetic field (B) induced electric po- 
larization or electric field (E) induced magnetic moment. 
The linear magnetoelectric response is characterized by 
the magnetoelectric tensor a and different physical 
mechanisms can contribute to it. One of these mecha- 
nisms has been recently the focus of considerable atten- 
tion d, H, 0, S [E 0] and is related to the concept of 
toroidic moment [8|. Materials with a uniform toroidic 
moment (ferrotoroidics) exhibit a non-vanishing magne- 
toelectric effect, and the ferrotorodic moment is an anti- 
symmetric component of a. These materials might have 
technological applications in the area of computer mem- 
ories. 

The observation of ferrotoroidic domains in LiCoP04 
confirmed the physical relevance of this concept and 
opened the possibility of observing toroidic driven mag- 
netoelectric effects in certain quantum magnets. How- 
ever, as it was pointed out recently Q, the toroidic mo- 
ment is multivalued in a periodic system. This problem 
is analogous to the case of the electric polarization. The 
root of this indeterminacy lies in the multivalued nature 
of the position operator for a periodic system: r = r -|- a, 
where a is a translation that leaves the system invariant 
0. The position operator appears explicitly in the ex- 
pressions for the multipolar moments of charge or current 
distributions and, thus, all of these moments are multi- 
valued in a periodic system. In particular, the toroidic 
moment appears in a multipolar expansion of a current 
density distribution and it couples linearly to V x B Q, 
and to E X B [l^ , which is easier to control. 

Since it is usually convenient to use periodic bound- 
ary conditions for modeling physical systems, it is neces- 
sary to know how to compute the ferrotoroidic moment, 
T, and the consequent magnetoelectric effect within this 
framework. Although r is multivalued for a periodic sys- 
tem, the change of r (At) between two different physi- 
cal states is well defined when moving along a particular 
path. In the case of the macroscopic electric polariza- 
tion, P, of an insulator, the change AP is computed by 



integrating the charge current through a given surface 
along an adiabatic path that connects the two different 
physical states d, [ll| . This leads to a natural relation 
between AP and the Berry phase associated to an adia- 
batic evolution in an enlarged parameter space pd] . [T3 | , 
and illustrates the observable character of the geometric 
phase. In the same way that the charge Berry phase 7^^ 
is a measure of the macroscopic electric polarization in 
band or Mott insulators [llj, the spin Berry phase 7" 
can be related to the difference between electrical polar- 
izations for spin up and down [Tsj . Here we will demon- 
strate that the change in the ferrotoroidic moment. At, 
can be related to a Berry phase tensor 7^^ , thus revealing 
the geometric character of the ferrotoroidic response. 

Ferrotoroidic moment. Assume a system of N inter- 
acting electrons of mass m and charge e enclosed in a 
d-dimensional box of linear dimension and volume 
il (/I = I,-'' jd). Particle j is labeled by the coordi- 
nates rj and spin Sj ^ \cf'^ {cr'^ are Pauli matrices with 
v ~ x,y, z). The Hamiltonian of the system is (ft = 1) 

N -rT2 N 

= ET^ + E^"'(l^'-^^l) + E^-t(r„A(i)) 
— ' Zm ^ — ' ^ — ' 

i—l i<Cj 2—1 

N 

+ 9so^'^H(l>iiu)(TiU {(t>p,u) -r-t X Hi, (1) 

1=1 

where Ilj = -I- A'^ is the canonical momentum (p^ = 
~idfj_) with gauge field A" = ^i^cr^e^ (e^ is a unit vector 
along the ^ direction). The flux is the twist in spin 
space at the boundary of the box: '^ai^i + L^e^j) = 
E^fe^^-'^la/J^Mr.). 

The transformation defined by U {(ppiu) = 
gi^f.^cr" r^/L,, g^^j^ ^^Yisx the eigenstates of 
satisfy periodic boundary conditions. The last term in 
(fT|) corresponds to the relativistic spin-orbit interac- 
tion, while Vint and Wgxt represent the interaction and 
"external" potentials, respectively. The latter depends 
parametrically on A(i) which is assumed to change 
adiabatically in time t from A(0) = to A(r) = 1. 

The second quantized form of is ~ 
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d^r TL^(y). with Hamiltonian density 



*t(r)IL<i,(r) 

gso*^(r)Ut 



^,)(tU (0^,).rxnvi/ (r),(2) 



written in terms of the fermionic spinor fields '4't(r) = 
('0|(r), '0x(r)), where V'l(r) (V'/3(r)) creates (annihilate) 
an electron with spin a{j3) = {Ti i} position r. 

The spin contribution to the ferrotoroidic moment 
is defined as the antisymmetric component of the tensor 
tpjy. The tpjy-density field is defined by the local condi- 
tion = — 25i,(r). For a finite system with 
open boundary conditions (surfaces), the following rela- 
tion holds as long as 7^y(r) vanishes outside the system 

= = ^ / d\ r, SAr) - ^ / d\ T,Ar), (3) 

written in terms of the spin density field 
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5.(r)-5^5^:^(r), 5^,,(r) 

ct/3 



V-KrX^V-Mr), (4) 



where ^^^d^r S^{y) = to have t^;^ independent of 
the origin of coordinates More specifically t,, = 

^Spi/^w^M'^; ■with £,,^1/ the Levi-Civita tensor. Like 
for the electric polarization [111, [l^, the third member 
of Eq. ^ is not well defined for a system with periodic 
boundary conditions because the position operator is not 
well defined [§]. However, the last member of Eq. ([3]) is 
still well defined, and starting from such expression for 
we will demonstrate that the change (At^(^) between 
two different physical states can be obtained from a set 
of Abelian Berry phases. 

Since the ground state (GS) of H'^ evolves continuously 
between A(0) = and A(T) = 1, the change is given by 



\\ f dVaA($^IV(r)|<I'^), (5) 
Jn 



where |$o) is the GS of with flux (^p^, and (At^^) = 

($^^^Vm'>I<I'o'^^) - {f^°htuAK'^°^)- By using adiabatic 
perturbation theory [111 . Ilq up to first order in the time 
derivative we obtain 



i^dV9,($o^lV(r)|$o^) = 



where we have used that 

ii'^inl^i') = S,n,m') and 



C.C., (6) 
(7) 



V(r)] = dtT^, = L^d4,^,^n\r), 

In the absence of the spin-orbit term, dfT,,,^ = 2J^jj{y) 
with J7^i^(r) the density of spin current [131 satisfying 



the continuity equation d^J^^,{r) + dtS^{r) = 0. By 

using the relations {<^^\dx^^\^o) = ('^'il [^a, H^] |$o)' 
and = we obtain 

from Eqs. © and dH) 
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»((aA<f>o^|a^,.$o'> - (a*..<I'o'|9A$o'))- Since 



is a bulk property (and we assumed that there is no 



with S(0 

level crossing for the GS as a function of </>^,y), its value 
should not depend on the boundary conditions when we 
take the thermodynamic limit. Thus [l3 | 



(Ai^,) 



with the line integral performed along the contour F 
of [0, 27r] X [0,1] in the plane ((?!>^y,A), and A{Q_ 



i($o I^C'^'o)- Assuming a path-independent gauge [ll| 



(AV) = ^ [7^.(1) -7^.(0)] 



with (A) 



/o 



(9) 



One can also introduce fluxes along space directions other 
than ^ and, as long as the gap does not close, average 
over these additional fluxes. The Berry phases 7^^, are 
anholonomies associated with the parallel transport of 
a vector state (GS) in the parameter space determined 
by the fluxes. Equation Q shows that the quantum of 
uncertainty^ related to the lack of history in the adiabatic 
evolution, is at least the inverse of the total transverse 
section L^/Cl. In a periodic system with primitive unit 
cell volume fio < i^, the quantum is larger and equal to 
L°/fJo, where L° is the length of the unit cell along the 
/i-direction [Tl| . 

The extension to lattice systems is straightforward 



rr' ,af3 



H.c. 



^fnt, (10) 



where cj^ creates an electron of spin a at lattice site r, 
The spin current on the bond (r,r + e^) is JT^^, = 



(ct 
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a'^g. The charge 
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Berry phase 7"^ is obtained when a'^ is replaced by 

drfj,Af^ = 0^(r^— rp)/L^, i.e., the flux associated with 
the usual vector potential of the electromagnetic field. 

Localization Indicators. There is a formal mathemati- 
cal connection between the geometric phases (7°, 7*) and 
the localization indicators (z£,2;£) @, [la, [3 that were 
introduced to discriminate between different phases such 
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as conductors and insulators |2C 



is defined by Eq. 



(fTTj) below, replacing lys^j.cr'^ with a scalar coupling (^^. 
For a one dimensional system along the x direction and 
spin quantization axis z, z£ corresponds to the compo- 
nent Zxz of the tensor defined by Eq. (fTTj) . The phase of 
the localization number z'l (with fc = c, s), is related to 
the Berry phase through the relation 7*^ = Im In z'l . 

Different indicators provide complementary informa- 
tion. They may display different convergence properties 
to the thermodynamic limit, and most importantly, some 
indicators are more akin to Monte Carlo methods while 
others are geared towards numerical renormalization 
group approaches. The localization parameters z^ and 
their related Berry phases j'^ have also been used as sharp 
topological indicators for establishing quantum phase di- 
agrams of interacting systems [HI, [l9|, [2l|, [l^, [H, [131 ■ 

By analogy to z^'", define the localization indicators 



(11) 

where ip^^ = ^n/L^ and Iv^o) is the GS [l4j . 

Considered as a continuous function of f^u^-, z^^^ plays 
the role of a characteristic function generating all mo- 
ments of the T^i/ tensor, and trivially satisfies 2;^,y[0] = 1, 
\z^4Vfj^i^]\ < 1: Zf^^[ipf^^] = z*^^[-Lp^y\. Note, however, 
that the operator T^jy is not a genuine operator in the 
Hilbert space bundle defined above, although its expo- 
nential is a legitimate one. Therefore, expectation values 
of arbitrary powers of T^,y have only meaning in terms 
of z^n (i.e., there is a quantum of uncertainty [IH). As- 
suming analyticity in the neighborhood of ip^i, = (i.e., 
the system has a gap in the thermodynamic limit), z^i, 
can be written in terms of cumulants Ck(T^u) 



exp 



f:^c.(r,.) 



Lfc=i 

with the end result that 

-1 



LP Imlnz^^ = {T^,^) + 0{lp ) , 



In I 



P = (TL) - {T^.f + 0{pI 



(12) 



(13) 
(14) 



so that the phase of z^j^ is related to the tensor T^j, while 
its modulus provides information on its quantum fluctu- 
ations (it vanishes in the thermodynamic limit for sys- 
tems with gapless excitations that have a non-zero spin 
current) . 

To illustrate an elementary application of Eq. 
assume that \^o) is the GS of a periodic system with 
the primitive unit cell of volume Oq, (such as that of 
LiCoP04 [^ but with spins pointing in the z direction) 
consisting of two spin up and two down in the direction 
z, and displaced in the direction x so that the contribu- 
tion to Tj-z from each unit cell, tj.^ = '^^j^iXjSj ^ 

((^^1^) — i^i//f^o)- Thus, the first member of Eq. 
should give N^t^i, where iV„ is the number of unit cells. 
This equation is invariant if any coordinate is re- 
placed by + nL^ with n integer. Equation pip gives 



z^^ = exp(27riiV„t^^/L^), and replacing into Eq. p3)) . 
(^Tti/) = Nyt^i, as expected. Note that the same result is 
obtained if t^,^ is replaced by t^,^ + nq, where n is inte- 
ger and q = Lfi/Nu is the quantum of uncertainty in i^^. 
However, for translational invariant systems, one can cal- 
culate z^j, in the effective one-dimensional problem with 
fixed total transversal wave vector, {Ky, Kz), for which 

Effective single-body schemes. One can implement pre- 
vious ideas in the framework of spin density (matrix) 
functional or Hartree-Fock theories. The Hamiltonian 
acting on the Kohn-Sham (or Hartree-Fock) orbitals, 
■05^j.(r), with lattice periodicity a 



, . - 

Ti— 





(£r'^"/3 + ^";3(r))V'fk(i") = enk<k(r) 



(15) 



defines the two-component spinor 5'„k(r) 



, ^ ) with 
the following generalized Bloch-spinor conditions 



*nk(r + a^e^) 
*nk(r) 



•i(k+:: 



'unk{r) (16) 



where M„k(r + a) = M„k(r) is the periodic part of the 
Bloch-spinor. In Eq. p^ Vap represents the effective 
one-body potential (which includes the external, Hartree, 
exchange and correlation components). To simplify the 
final expression of {t^u), it is convenient to express the 
spinor it„k in the basis of eigenstates of a'^ , 



, with a^Unu 



(17) 



and the periodic part of the (non-interacting) many- 
particle wavefunction along each of the spin directions 
is the product of Slater determinants: $J^(ri, • • • , tat) = 
Det(A^+(k))Det(A''-(k)), where Ar/(k) = <k^(rj) with 
e = ±, 1 < i < Ng, Siud 1 < j < N+ (1 + N+ <j <N) 
for 6* = + (6* = -). N ^ N+ + is the number of 
particles in the unit cell. Then, using Eq. ([9]) 



(V(A)) 



(27r)3 



/ rf'fc^«ki 5fc.<k),(i8) 

9=± n=l -^^2 



with the momentum integral evaluated over the Brillouin 
zone (BZ) corresponding to the periodicity a. As men- 
tioned above, what has physical meaning is the change 
in the tensor t^^. Therefore, At^i, = i^iy(l) — t^,y{0). 
Here, as in the case of charge and spin macroscopic po- 
larizations [m , the price paid for considering a two-point 
formula (i.e., forgetting about the A-dependent path) is 
the appearance of a quantum of uncertainty. 

Geometry of magnetoelectric response. Finally, to un- 
derstand the geometric content of the magnetoelectric re- 
sponse, we need to understand the Riemannian structure 
of our Hilbert space bundle. Consider a set of normalized 
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states {l^'oCv))}: where represents 3d real numbers 
(/j^y. Let's assume that this manifold of quantum states 
is generated by the action of the group of transformations 
l^oif)) = e'¥'-T|v[/o(0)), where ^-T = E^. ^^m-T^a.- 
The expectation value of the generalized twist operator 
e'V T measures the character of the spectrum of low- 
energy T^i^-excitations. Can we find a measure of the 
distance between two of these quantum states? 

Provost and Vallee addressed the problem of estab- 
lishing a Riemannian structure on an arbitrary differen- 
tiable manifold of quantum states psj . Berry, later on, 
extended their work by introducing a geometric tensor 
[2^ whose physical relevance in the context of electric 
charge polarization was addressed in Refs. [H, [l^l and 
also for spin polarization in Ref . [3 ■ We now introduce 
a general quantum geometric tensor 

G^,;M = (5^„.*o|Po[<p]|a^„,*o), (19) 

whose real and imaginary parts G^^-ap[v] ~ 9fj,iy;ap[v] + 
i ^fj,}^;af3['p] can be written 

{ (20) 

with Pa[ip] = 1 - |*o(¥'))(*o(¥')L and "/^^if] = 
*(^o|(?ip,,„^o) the Berry connection [ll|. The real part 
of G i^i^-a/ii'P] is a symmetric and positive definite ten- 



sor representing a generalization of the metric intro- 
duced in Ref. [2^. Moreover, it is interesting to remark 
that the infinitesimal distance is related to the quantum 
fluctuations of the T^^, tensor, i.e., g^i^.a^i^ — > 0) = 
{T^^Tap) — {T^^){Tap), with expectation values evalu- 
ated over l^'o(O)) = |^'o)- In a sense, the metric struc- 
ture on the manifold is fixed by the quantum fluctua- 
tions which determine the modulus of z^^ in the ther- 
modynamic limit. On the other hand, the antisymmet- 
ric tensor n^^.afilv] = Ini(5^^„*o|l9v„3^o) {^tiy;apW\ = 

^^a/3■,^ll'Vp\) plays the role of a curvature, and is a quan- 
tity connected to the non-dissipative part of the spin con- 
ductance in adiabatic transport. 

In summary, we have introduced a formalism that not 
only leads to a geometric understanding of a spin compo- 
nent of the magnetoelectric response but also provides a 
computational method to study changes in ferrotoroidic 
moment. At, when the system under consideration has 
periodic boundary conditions (no surface). Forgetting 
about the history of the evolution leading to At amounts 
to the appearance of a quantum of uncertainty of magni- 

tude ^ Our many-body formalism generalizes, giv- 
ing an operational physical interpretation to, the concept 
of spin Berry phase introduced in [isj . and should be 
particularly useful within the framework of spin density 
(matrix) functional or Hartree-Fock theories. 

LANL is supported by US DOE under Contract No. 
W-7405-ENG-36. 
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